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INTRODUCTION 

Analysis  of  the  delay  experienced  by  a  side-road  driver  attempt¬ 
ing  to  merge  with,  or  cross,  traffic  on  a  main  road  has  been  conducted  by 
many  authors*  see  Bisbee  and  Oliver  [  2  ],  Buckley  and  Blunden  [3],  Evans, 
Herman,  and  Weiss  [4],  Garwood  [5],  Gaver  [6],  Hawkes[  8],  Jewell [9],  and 
others.  Postulating  various  behavioral  characteristics  for  side-road 
drivers,  the  above  authors  have  presented  formulas  describing  the  long- 
run  delays  of  drivers  attempting  to  merge  or  cross,  and  other  related 
measures  such  as  side-road  queue  magnitudes. 

The  purpose  of  the  present  paper  is  to  study  one  such  model 
numerically,  with  special  emphasis  upon  time -dependent  behavior.  For 

particular  examples  we  shall  exhibit  the  manner  in  which  expected  delay 
either  converges  to  the  long-run  limit,  if  a  finite  limit  exists,  or 
grows  if  there  is  saturation.  Without  such  information  we  are  forced  to 
rely  upon  the  "steady-state"  or  long-run  results  of  queueing  theory. 
Steady- state  values  may  be  reached  rather  slowly  under  some  conditions, 
e.g.  during  rush  hours*  it  is  perhaps  useful  to  have  some  idea  of  the 
adequacy  of  the  long-run  predictions  at  various  times  after  the  start  of 
the  rush  hour.  One  factor  that  emerges  as  important  in  this  study  is  the 

effect  of  mixing  "fast"  with  "slow"  diivers  in  various  proportions  on  the 
If  I 
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aide-road,  "feat”  drivera  being  thoae  willing  to  accept  relatively  ahort 
main-road  gapa  for  merge  or  crossing,  "alow  drivers  (e.g.  possibly 
buses  and  trucks)  requiring  somewhat  more  time. 

Finally,  our  numerical  procedure  allows  the  study  of  natural 
approximations  to  the  queueing  processes  that  arise  in  this,  and  other, 
contexts.  We  give  an  example  of  one  such  approximation  that  is  related 
to  diffusion  theory  and  the  "heavy  traffic"  approximations  of  J.F.C. 
Kingman.  It  seems  sensible  to  attempt  to  find  plausible  and  comprehen¬ 
sible  approximations  for  the  queueing  processes  arising  in  traffic,  as  in 
other  areas  of  applied  probability.  Numerical  methods  provide  one  means 
of  discovering  and  evaluating  such  approximations,  and  of  uncovering 
important  differences  in  the  various  models  that  have  been  proposed. 

THE  MODEL 

Our  basic  model  haa  been  introduced  in  an  earlier  paper  [6], 

We  shall  assume  that  main  road  traffic  is  a  stationary  Poisson  process 
with  rate  v  and  that  main  road  vehicles  can  be  considered  to  be  points. 
Side-road  or  merging  traffic  is  likewise  Poisson  with  rate  X.  There 
are  two  types  of  side-road  drivers:  Type  I  will  be  called,  for  want  of  a 
better  term  "fast"  and  Type  2  will  be  called  "slow",  although  other 
adjectives  may  be  more  meaningful)  the  restriction  to  two  types  is 
easily  lifted.  The  probability  that  fast  and  slow  drivers  appear  at  the 
merger  of  the  roads  in  any  order  is  that  of  a  sequence  of  Bernoulli  trials: 
the  probability  of  a  Type  1  being  p^  and  of  a  Type  2  being  p2<  The 
reasons  for  distinguishing  between  the  two  driver  types  is  that  each  type 
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has  a  distinct  gap  acceptance  mechanism;  typically  a  fast  driver  will  accept 
or  merge  into  a  shorter  main  road  gap  than  will  a  slew  driver,  so,  condi¬ 
tional  upon  the  driver  type,  his  gap  acceptance  probability  function  will 
be  given  and  he  will  behave  in  accordance  with  it.  Putting  this  another 
way,  a  side  road  driver  experiences  interrupted  service:  he  compares  his 
critical  gap  to  each  new  gap  in  ma^a  road  traffic.  In  our  present  model 
he  selects  a  new  critical  gap  from  his  characteristic  gap  acceptance 
probability  distribution  for  comparison  to  each  new  main  road  gap.  When 
for  the  first  time  his  critical  gap  is  less  than  the  main  road  gap  he 

makes  entry  into  the  main  road  stream.  Furthermore  —  and  this  is  a 

\ 

distinctive  feature  of  this  model  —  the  side  road  driver  indeed  uses  the 
entirety  of  his  critical  gap  while  making  entry.  In  actual  fact,  drivers 
perhaps  "play  safe",  and  wait  for  a  main- road  gap  that  is  much  longer  than 
necessary  for  entry.  However,  if  main  road  traffic  is  reasonably  heavy 
and  gaps  typically  short,  drivers  may  well  behave  approximately  as  postulated 
here.  If  so,  and  if  main-road  traffic  is  Poisson,  then  instants  of  entry 
from  the  side  road  are  regeneration  points,  and  the  side-road  queue  is  an 
ordinary  single-server  system  with  independent  completion  times  in  place 
of  service  times.  For  simplicity  we  shall  consider  the  service  (completion) 
time  to  be  made  up  entirely  of  the  wait  at  the  merge  point  (head  of  the 
queue)  for  a  suitable  gap  to  appear,  and  neglect  the  move-up  time.  This 
would  seem  reasonable  under  the  conditions  for  which  the  present  model 


is  applicable. 
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GAP  ACCEPTANCE  AND  COMPLETION  TIMES 

Our  aim  is  to  give  numerical  results  for  the  model,  so  we  have 
assumed  a  specific  gap  acceptance  probability  function  for  both  "fast"  and 
"slow"  drivers.  The  probability  density  of  the  critical  gap  for  a  driver 
of  "type  i  (i=l(fast),  =$(slov))  is  gamma: 


t  k,  k.-l  .  t  k.  k 

fi(x>  v  V  =  fry  expt-  ^  • 


(i) 


Since  the  expected  value  of  the  critical  gap  is  m^,  and  its  variance  is 
2  -1 

mi  ^i'  cr^^c£*^-  8aPa  a  particular  side-road  driver  may  be  forced 

to  cluster  as  tightly  as  desired  around  any  value  for  m^  by  simply 
increasing  k^.  We  shall  shortly  illustrate  the  effect  of  changing  the 
variance  of  the  critical  gap  while  its  mean  is  held  fixed. 

The  gap  acceptance  procedure  is  as  frllows.  When  a  new  driver 
moves  into  merge  position,  he  is  of  type  i  with  probability  p^(i=l,2).  He 
immediately  selects  a  critical  gap,  G^,  from  density  f^.  If  there  are  no 
main-road  cars  to  arrive  at  the  merger  within  time  G^  he  enters,  completing 
entry  at  time  0^  after  he  initially  reached  merge  position.  Suppose,  however, 
that  a  main-road  car  will  reach  the  entry  within  G^.  The  side-road  driver 
waits  until  the  instant  it  passes  the  merge  point,  and  then  selects  another, 
independent,  gap,  G2,  from  (l),  and  behaves  as  he  did  with  G^,  and  so 
forth  for  G^,  etc.  The  probability  that  a  side-road  driver  enters  during 
any  critical  gap  is  the  probability  that  the  latter  contains  no  main-road 
car,  i.e. 

r°  vm  -k 

cp±( v)  =  /  e-vx  f±(x,  k±,  mi)dx  =  (l  +  -~) 


(2) 
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and  the  Laplace  transform  of  the  density  function  of  a  critical  gap, 
given  that  entry  is  made,  is 


oo 


Cpi(\M-s) 


1  r  -ax  -vx  „  ,  ,  V1"8 

5^7  le  e  fi(x>  V  "i)ix  = 


(3) 


-sO' 


If  O'  denotes  the  above  conditional  critical  gap,  then  E[e  ]  is  given 
by  (3).  The  Laplace  transform  of  the  density  function  of  a  main  road 
gap,  given  that  it  is  shorter  than  the  critical  gap  and  hence  no  entry 
is  made,  is,  putting  F^x)  for  the  distribution  function  associated  with 
(1), 

i  r  -sx  f  ,  n  -vx  i“<Pj(8+v)  v 

- - t — r  •  /  e  [l-F.(x)Je  v  dx  =  — - r-t - .  (4) 

1“<P±C  v)  ^  1  iWJ  1-^Tv)  v+s 

-bM 

If  M  denotes  such  a  conditional  main-road  gap,  then  E[e  ]  is  given  by 
(4).  Finally 

c  =  m1  +  m2+...  +  m^+o*  (5) 

where  N  denotes  the  number  of  too-short  main  road  gaps;  clearly 

?|N=n|Type  i|  =  [l-cpi( v)  Jn  <p±( v)  .  (n=0,  1,  2,  ...)  (6) 

It  then  follows  easily  from  the  apparent  independence  that 

oo  n 


E[e  13  =  |E[e"sM]|  E[e'sG  ] [l-cp^^C  v)  ]ncp±(  v) 

n=0 


(a+vjcp^s+v) 

s+vcp^s+v) 


(7) 


The  Laplace  transform  of  the  completion  time  of  a  random 


driver  attempting  to  merge  is 

2 


_sC  (3+v)<Pi(8+V) 

C  ^  pi  ,3+vcpi(s+v)  ' 


(8) 
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I 


and 


E[CS] 


l-<p±(  v)h  vq>j(  v) 

2  2,  > 

V  <P±(  v) 


(9) 


q>^(v)  denotes  the  first  derivative  of  <p^  at  v.  Nov  the  long-run  ex¬ 
pected  waiting  time  for  a  side-road  driver  to  enter  can  be  obtained  under 
our  assumptions  by  using  C  as  the  service  time  in  the  classical  formula: 


^  ■  51^7  IB?  '  (w) 

provided  p  =  XE[C]  <  lj  here  W(t)  denotes  waiting  time  at  an  irstant  t 

time  units  from  process  Initiation. 

The  expression  (10)  is  limited  to  assessing  expected  waits  "a 

long  time"  after  the  process  begins,  and  then  only  when  p  <  1.  Thus  it  is 

of  no  direct  use  for  predicting  delays  at  any  specific  short  time  after 

some  initial  instant  —  at  which  conditions,  e.g.  v  or  X,  changed,  as  at 

the  onset  of  rush  hour  --or  when  the  merging  capacity  is  saturated,  i.e. 

when  p  >  1.  In  order  to  exhibit  the  manner  in  which  expected  waiting  time 

depends  upon  elapsed  time,  t,  and  also  upon  the  initial  state  cf  the 

process,  one  way  to  proceed  is  to  numerically  invert  the  Laplace  transform 

of  the  expected  waiting  time.  Specifically,  letting  N(t)  denote  the 

number  of  cars  queued  at  the  side  road  at  time  t,  we  invert 

oo 

p( s )  =  £  e"8t  E[W(t)|N(0)  =  i]dt,  (ll) 

a  formula  for  which  may  be  derived  from  general  queueing  theory.  Hiis 
formula  is  presented  in  [7],  where  a  method  of  numerical  transform  inversion 
is  also  developed.  The  methods  of  [7]  have  been  applied  to  several  specific 
cases  to  obtain  the  numbers  in  the  following  table. 


\ 
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EXPECTED  DELAY  AT  SIDE  ROAD 
(initially,  no  queue) 

Px  =  0.5,  p2  =  0.5 

m1  =  5#  m2  =  15  (secs) 


Case 

\f(sec8) 

k 

60  (secs)  12.0 

1  =  k2 

180 

=  2 

240 

600 

900 

Long-Run 

do) 

I 

X  *  0. 05 

V  =  0.10 

E(S)  =  13.12 
p  =  O.656 

(14.6 

27.7 

19.5 

32.6 

22.2 

35.3 

24.0 

37.1 

27.7 

40.8 

28.3) 

41.4 

(28.7) 

41.8 

II 

X  =  0.05 

v  =  0. 15 
e(s)  =  14.69 
p  =  0.735 

(18.4 

32.1 

25.8 

4o.5 

30.5 

45.2 

33.9 

48.6 

43.1 

57.8 

45.9) 

60.6 

(49.1) 

63.8 

III 

X  =  0. 075 
v  =  0.10 
E(S)  =  13.12 
p  *  0.984 

(25.7 

38.8 

39.0 

52.1 

49.2 

62.3 

57.7 

70.8 

94.1 

107.2 

115.8) 

128.9 

(946.9) 

960.0 

IV 

X  =  0. 075 

v  =  0.15 
E(S)  =  14.69 

p  =  1.102 

(31.9 

U6.6 

50.4 

65.1 

65.6 

80.3 

78.9 

93.6 

143.2 

157.9 

187.9) 

202.6 

Not 

Applicable 

P1  = 

0.25,  p2  * 

0.75  (otherwise  same 

as  above) 

V 

x  =  0.05 

v  =  0.10 
E(S)  =  16.68 
p  =  0. 844 

(22.1 

39-0 

32.0 

48.9 

38.8 

55-7 

44.0 

60.9 

61.9 

78.8 

69.7) 

86.6 

(90.6) 

107.5 

P1  = 

0.75,  P2  * 

0.25  (otherwise  same 

as  above) 

VI 

X  =  0.05 
v  =  0. 10 

E(S)  =  9.38 
p  =  0. 469 

( 7.8 
17.2 

9.4 

18.8 

10.0 

19.4 

10.2 

19.6 

10.5 

19.9 

10.5) 

19.9 

(10.5) 

19.9 

a. 


EXPECTED  DELAY  AT  SIDE  ROAD( Continued) 
(initially!  no  queue) 


Case\t(seca) 


P,  =  0.5,  PP  -  0.5 


(otherwise  same  as  above) 

600  900 


Long-Run 


VII 


x  *  0.05 

v  =  0. 10 
E(S)  0  17.15 
p  -  O.858 

(23.3 

Vo.  5 

34.2 

51.4 

41.9 

59.1 

48.0 

65.2 

69.5 

06.7 

79.5) 

97.7 

(112.7) 

129.9 

X  =  0.075 

v  =  0.10 

I(32  “1 2&15 

(40.4 

57.6 

66.7 

83.9 

89.8 

107.0 

112.2 

129.4 

226.6 

143.8 

316.6) 

Not 

Applicable 

VIII 


NOTE:  Numbers  in  parentheses  indicate  expected  wait  until  head  of  line  is 
reached j  othsr  numbers  Indicate  total  time  to  make  entry. 

Dimensions  of  X  and  v  are  [secs"lj,  and  of  E(S)  are  [secs]. 


r 


‘“f 


■i H’-'  mif 
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Although  Table  1  was  computed  for  specific  parameter  values,  It 

helps  to  provide  an  understanding  of  the  approximation  afforded  by  the  long 

run  formula,  (10),  when  the  latter  is  applicable.  For  comparison,  the  long 

run  values  of  expected  waiting  time  corresponding  to  the  various  oases 

appear  in  the  right-most  column  of  the  table.  For  example, 

Case 

1)  Comparlng/I  to  Case  VII  suggests  the  important  effect  of  the 
"shape"  or  "sharpness"  of  the  gap  acceptance  probability  function,  as 
measured  by  the  parameter  k,  both  upon  the  long-run  expected  wait  and 
upon  the  rate  at  which  the  long-run  value ;  is  approached.  Although  the 
expected  critical  gaps  for  both  driver  types  are  the  same  in  both  cases,  in 
Case  VII  the  probability  that  a  side-road  driver  has  a  small  critical  gap 
is  much  smaller  than  in  Case  I.  In  the  present  model  this  is  because  the 
variance  of  the  critical  gaps  for  Case  VII  is  one-third  that  in  Case  I. 

The  result  is  that  entry  is  more  difficult,  and  the  traffic  intensity  is 
larger,  in  Case  VII,  and  the  long-run  value  is  much  less  rapidly 
approached. 

2)  Comparing  Cases  I,  V,  and  VI,  the  effect  of  driver  type  mix  is 
exhibited.  Clearly  a  mixture  of  25  percent  slow-75  percent  fast  results 
in  a  drastically  reduced  wait,  especially  when  compared  to  75  percent 
slow-25  percent  fast.  Moreover,  approach  to  the  long-run  value  is 
correspondingly  effected. 

3)  Finally,  Cases  III  and  VIII  indicate  the  effect  of  the  gap 
acceptance  probability  shape  upon  wait  when  the  merge  is  close  to 
saturation  (former  case)  or  over  saturation  (latter). 
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APPROXIMATIONS 

Examination  of  expressions  (5),  (6),  and  (7)  reveals  that 
la  a  geometrically  compounded  random  variable.  If  E(C^)  becomes  large, 
either  because  main-road  traffic  density  la  high,  or  the  probability  of  an 
entry  from  the  side-road,  ( v) ,  Is  small,  then  the  distribution  of 

la  approximately  exponential; 


P  °1  <  *]•  ^  1  •  -x/B [Cj  J  j-  (x  >  0). 


(12) 


Tendency  to  the  exponential  In  distribution  can  be  shown  formally  by 
means  of  a  continuity  theorem  for  Laplace  transforms. 

The  fact  that  the  service-time  distribution  tends  to  the 
exponential  suggests  that  an  approximating  process  can  be  constructed, 
having  exponential  service  times  in  place  of  the  true  service  times  of  (8), 
that  mimics  closely  both  the  transient  and  long-run  behavior  of  the  original 
process.  It  Is  this,  and  allied,  approximations  that  ve  shall  study 
numerically  here. 

Suppose  jw(t),  t  >  oj-  represents  the  original,  or  "true",  delay 
process  at  the  merge  point)  this  process  has  arrival  rate  X  and  completion 
times  described  by  (8).  In  the  light  of  (12),  a  first  candidate  for  an 

approximate  process  is  where  the  latter  has  arrival  rate  X  and 

mixed  exponential  completion  times: 

-oC  pi  pp 

*  =  1+sElcJ  +  i+se(c2J  * 

A  second  candidate  is  |w8(t)^-,  differing  from  |  W  ( t)|  only  In  replacing 
(13)  by  the  single  exponential 


11. 


-•c 

■  (•  "I 


1 _ 

PjMCj  ♦  PgElGjj 


(14) 


'Hie  degree  of  approximation  thus  obtained  la  Illustrated  In  the 


next  table  for  two  cases. 


IHble  2 

COMPARISON  Or  APPROXIMATIONS 
FOR  EXPECTED  DELAY 
(no  Initial  queue) 

A  »  O.C5,  v  -  O.lOj  m^  •=  5,  ®2  -  15 

-  *2  -  2j  E[C1 )  =  5.625,  E[C2 J  -  20.625 


I:  Pl  -  0.75,  P2  -  0.25*  E[C)  >  9.375 
P  *  0.469 

CaseW sees)  60  120  l6o  240  6 00  900 


E(W(t)1 

7.6 

9-4 

10.0 

10.2 

10.5 

10.5 

E(WB(t)] 

8.6 

10.6 

11.1* 

11.8 

12.2  • 

12.2 

*[«.(*)] 

6.9 

7.6 

8.1 

8.2 

8.3 

8.3 

II: 

P1  * 

0.25,  p2 

*  0. 75)  E[C J  *  16.875 

P  - 

0.844 

E[W(t)] 

22.1 

32.0 

38.8 

44.0 

61.9 

69.7 

E(W  (t) J 

UA 

23.3 

34.2 

41.7 

47.6 

68.1 

77.3 

E  [W#  (  T  )  J 

22.3 

32.1 

39.0 

44.2 

62.2 

70.1 

Apparently,  the  approximation  is  only  fair  in  I  and  respectable 
In  II,  with  the  single  exponential  approximation  superior  to  the  mixed 
exponential  in  II. 
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It  is  noticeable  in  I  that  neither  approximation  is  approaching 
the  limiting  value  for  the  true  process.  Unless  main-road  traffic  is  very 
dense  (v  ->  oo )  the  exponential  (12)  will  not  be  precise,  so  ve  cannot  in 
general  expect  the  approximating  processes  to  possess  the  same  long-run 


limits  as  does 


The  latter  limits  depend  upon  higher  moments  than 


the  first)  parenthetically  the  close  agreement  of  E[W  ( t) )  and  E[W(t)]  in 
II  above  can  be  explained  by  the  fact  that  the  second  moments  of  the 
respective  completion  times  are  very  close  numerically. 

The  source  of  an  Improved  approximation  is,  thus,  to  force  the 
approximating  process  to  have  the  same  long-run  mean  values  as 


as  veil  as  the  same  traffic  intensity.  Consider  (Wd(t))-,  having  arrival 


’}■ 


fw(t)} 


rate  and  exponential  completion  time  Cd  such  that 


XE[C  ] 


and 


E[W] 


2^1  ffr1  -  E(0dJ  -  E‘wd)  • 


(15) 


(16) 


Thus  ve  put 


(17) 


El°dJ  c  filci  1 

Another  vey  of  arriving  at  the  approximation  |wj(t)^-  is  to  argue 
in  a  manner  reminiscent  of  that  used  to  obtain  the  diffusion  approximation 
to  discontinuous  processes)  see  Bailey  (lj  or  Khintchine  [10].  Consider 
the  increments 


£W(t)  =  W(t+At)  -  W(t) 

£Wd(t)  =  Wd(UAt)  -  Vd(t)  . 


(19) 
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We  expect  the  development  of  W(t)  to  be  probebillatlcally  •ltnilar  to  that 
of  W^(  t)  If  the  distributions  of  the  Independent  increments  are  similar. 

To  bring  about  this  similarity,  equate  the  infinitesimal  means  and  variances 
of  the  tvo  processes)  neglecting  the  effects  of  the  boundary  at  aero  ve  have 


11m 

At 


11m 

AtrM) 


XE  [C]  -  1 
Xd  E[Cd]  -1 


Ha  4r  |  •  XE[C2) 

At-*0  ° 


llm  TrV«r[AH(t)J 

At-XO 


K  Etc/)  . 


(20) 

(21) 

(22) 

(23) 


Now  equating  traffic  intensities  in  the  tvo  processes  is  precisely  equiva¬ 
lent  to  equating  the  Infinitesimal  means  or  "drifts"  of  (20)  and  (21). 

Then  if  we  specialise  Cd  to  have  the  exponential  distribution,  and  equate 
(22)  and  (23)  the  expressions  (17)  and  (l8)  fall  out. 

The  qualities  of  the  tvo  approximations  suggested  are  illustrated 
in  the  following  table. 

Table  3 

COMPARISON  OF  APPROXIMATIONS  FOR  EXPECTED  DELAY 
(No  initial  queue) 

pl  «  P2  -  0.5  X  -  0.015,  v  -  0.25)  E[C ]  -  39.051 

=  5,  nt,  »  15  Xd  -  0.01016,  E (Cd )  -  57.665 


Cas^r(secs) 

60 

130 

180 

240 

600 

.900 

1200 

E[W(t)J 

22.45 

34.81 

43.09 

49.17 

67.23 

73.31 

76.54 

E[W8(t)J 

20.41 

29.68 

35.41 

39.36 

49.80 

52.62 

53.89 

E[Wd(t)] 

23.44 

35.91 

44.15 

50.16 

67.84 

73.73 

76.83 

14. 


Although  the  approximation  afforded  by 


fw.<4 


ia  inadequate, 


that  given  by  jwcJ(t)|  would  probably  be  entirely  aatiafactor  for 
applications. 

At  present  it  is  not  known  how  well  a  direct  diffusion 


approximation  to 


matlon 


fwd(t)}. 


fw(t)} 


compares  in  quality  to  the  exponential  approxi- 


The  diffusion  approximation  can  be  obtained,  in  principle. 


by  solving  the  forward  partial  differential  equation  of  diffusion  theory 
for  the  density  of  waiting  time,  and  then  Integrating  to  find  the  mean. 
Neither  is  it  yet  known  how  closely  higher  moments,  e.g.  the  variance, 


h(t>} 


conform  to  corresponding  moments  of 


fw(t)}, 


nor  how  well  the 


probabilities  of  waits  of  various  durations  in  the  two  processes  agree. 
J.F.C.  Kingman's  heavy  traffic  theory  [11],  [12]  does  suggest  that,  in  case 
side -road  arrival  rate,  X,  Increases  so  that  traffic  intensity  approaches 
unity  from  below,  the  diffusion  approximation  becomes  increasingly 


accurate.  An  increase  in  main-road  density,  v,  has  a  similar  effect. 


15. 
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